Quantum Phase Transitions and Matrix Product States in Spin Ladders 

M. AsoudehQ V. Karimipour^^Jil and A. Sadrolashraf 



Department of Physics, Sharif University of Technology, P.O. Box 1 1365-9161, Tehran, Iran 

(Dated: February 1, 2008) 

Abstract 

We investigate quantum phase transitions in ladders of spin 1/2 particles by engineering suitable matrix 
product states for these ladders. We take into account both discrete and continuous symmetries and provide 
general classes of such models. We also study the behavior of entanglement between different neighboring 
sites near the transition point and show that quantum phase transitions in these systems are accompanied by 
divergences in derivatives of entanglement. 

PACS numbers: 75.10.Jm, 73.43.Nq, 75.25.+Z 



^ Corresponding Author 

*Electronic address: asoudeh @ mehr. s harif .edul 
^Electronic address: vahid@ sharif. eduj 
^Electronic address: Isadr ' of nobles @ sharif .edul 

1 



I. INTRODUCTION 



The basic paradigm of many body physics is to use analytical and numerical tools to investigate 
the low-lying states and in particular the ground state properties of a system governed by a given 
Hamiltonian. At very low temperatures, when thermal fluctuations are dominated by quantum 
fluctuations, quantum phase transitions can occur due to the change of character of the ground 
state [HI]. What is exactly meant by "character" has been investigated in numerous works [E!], [jsj], 
[0], especially in recent years after the discovery of exact measures [Q] of entanglement(purely 
quantum mechanical correlations). For example, it has been a matter of debate whether a quantum 
phase transition is always accompanied by a divergence of some property in entanglement of the 
ground state wave function. 

Unfortunately, except for a few exactly solvable examples, the task of finding the exact ground 
state of a given Hamiltonian is notoriously difficult. As always in dealing with difficult problems, 
one way round the difficulty is to investigate the inverse problem, that is to start from states with 
pre-determined properties and investigate quantum phase transitions which occur by smoothly 
changing some continuous parameters of these states. The suitable formalism for following this 
path is the matrix product formalism [l6|], [ItI], [B, Is i which in recent years has been followed in 
constructing various models of interacting spins [|ld . 1151] . In this paper we want to 

take one step in this direction and in particular, we want to construct as concrete models, ladders 
of spin one-half particles and see what happens to the entanglement between various spins when 
the system undergoes a phase transition. We construct general class of models having a number 
of discrete and continuous symmetries. These are symmetry under the exchange of legs of the 
ladder, symmetry under spin flip, and symmetry under parity (left-right reflection of the ladder) 
in addition to a continuous symmetry, namely rotation of spins around the z axis or all three axes 
(full rotational symmetry). 

In these models we calculate the entanglement of one rung with the rest of the ladder as measured 
by the von Neumann entropy of the state of the rung and the entanglement of the two spins of a 
rung with each other as be measured by the concurrence of the same state. 

We will see that in all these models quantum phase transitions occur at a critical point of the 
coupling constant, and this point is where a rung of the ladder becomes completely disentangled 
from the rest of the ladder and the two spins of the rung become fully entangled with each other. 
The derivatives of these two types of entanglement are also divergent. 

We should stress here that we are using the term phase transition in a wider sense than usual Isj], 
that is we call any discontinuity in an observable quantity (i.e. a two point correlation function), a 
phase transition. 

The structure of this paper is as follows. In section [III we review the formalism of Matrix 
Product States (MPS) with emphasis on the symmetry properties of such states. In section Hill 
we specify the MPS construction to ladders of spin 1/2 particles and set the general ground for 
construction of concrete models. In this same section we construct multi -parameter families of 
models which have specific symmetries. In sectionHV] we study in detail the properties of the con- 
structed states and in particular calculate the exact correlation functions of spins on the rungs. We 
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also investigate the connection between quantum phase transitions and divergence of entanglement 
properties in these models. Finally we derive in the appendix, the Hamiltonian which governs the 
interaction of these models for which the state we have constructed is the exact ground state. 



II. MATRIX PRODUCT STATES 



First let us review the basics of matrix product states [0, S]. Consider a homogeneous ring of 
N sites, where each site describes a (i— level state. The Hilbert space of each site is spanned by the 
basis vectors \i), i = 0, ■ ■ ■ d — 1. A state 



iivPl; "^2, ■ ■ ■ , In) 



(1) 



is called a matrix product state if there exists D dimensional complex matrices Ai E 
■ ■ - d — 1 such that 



■'DxD 



1 



where Z is a normalization constant given by 

Z = tr{E^) 

and 

d-l 



(2) 



(3) 



(4) 



i=0 



Here we are restricting ourselves to translationally invariant states, by taking the matrices to be 
site-independent. For open boundary conditions, any state can be MPS -represented, provided that 
we allow site-dependent matrices Af^\ where k denotes the position of the site EQ. The 
MPS representation ^ is not unique and a transformation such as UAiU~^ leaves the 

state invariant. In view of this we can find the conditions on the matrices which impose discrete 
symmetries on the state. The state is reflection symmetric if there exists a matrix U such that 
Af = UAiU~^ where is the transpose of A and time-reversal invariant if there exist a matrix 
F such that 4* = VAiV-^. 

Let O be any local operator acting on a single site. Then we can obtain the one-point function on 
site k of the chain (\&|0(A;)|\&) as follows: 



mo{km 



tr{E- 



^-^EqE^-^) 



tr{E 



where 



d-l 



Eo:= J2mj)A*®Aj. 

i,j=0 

In the thermodynamic limit oo, equation ([5]) gives 

( ^max I -^O I ^max ) 



(^|0|^) 



A. 



(5) 



(6) 



(7) 
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where we have used the translation invariance of the model and Xmax is the eigenvalue of E with 
the largest absolute value and \Xmax) and {Xmax\ are the right and left eigenvectors corresponding 
to this eigenvalue, normalized such that {Xmax\Xmax) = 1- Here we are assuming that the largest 
eigenvalue of is non-degenerate. 

The n-point functions can be obtained in a similar way. For example, the two-point function 
(\[^|0(A;)0(/)|\&) can be obtained as 

(,io(.,o(oi.) ^ (8) 

where Eo{k) := E'^^^EqE^'^. Note that this is a formal notation which allows us to write the 
n-point functions in a uniform way, it does not require that E is an invertible matrix. In the 
thermodynamic limit the two point function turns out to be 

(^|0(l)0(r)|^) = -^Y.^'i''^^"---\^o\\){MEo\)^max). (9) 



max 



For large distances r ^ 1, this formula reduces to 



> r-2 



(^|0(l)0(r)|^) - {^\0\^f = ^{X^^^\Eo\Xi){X,\Eo\X^ax), (10) 



'max 



where Ai is the second largest eigenvalue of E for which the matrix element {Xi\Eo\Xmax) is non- 
vanishing and we have assumed that the eigenvectors ofE' have been normalized, i.e. (Ai|Aj) = 5ij. 
Thus the correlation length is given by 

I y-. ^Tnax 
Al 

Any level crossing in the largest eigenvalue of the matrix E signals a possible quantum phase 
transition. Also, due to (fTTl) . any level crossing in the second largest eigenvalue of E implies 
the correlation length of the system has undergone a discontinuous change. Here we are using 
a broader definition of quantum phase transition, that is we call any non-analytical behavior of a 
macroscopic property, a quantum phase transition Of course in the models that we construct 
we observe a more direct change of observable physical properties, namely in one regime we have 
correlations between spin operators on different sites and in the other we have no such correlation. 



A. Symmetries 

Consider now a local symmetry operator R acting on a site as R\i) = Rji\j) where summation 
convention is being used. Ris a d dimensional unitary representation of the symmetry. A global 
symmetry operator TZ := i?®^ will then change this state to another matrix product state 

vl/,,,,..,,-.M/':=tr(A: A:^---A:J, (12) 
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where 

A, ■= R,,Aj. (13) 

A sufficient but not necessary condition for the state |\E') to be invariant under this symmetry is that 
there exist an operator U (R) such that 

R,jAj = U{R)AiU~\R). (14) 

Thus R and U{R) are two unitary representations of the symmetry, respectively of dimensions d 
and D. In case that i? is a continuous symmetry with generators Ta, equation (fT4)) . leads to 

{Ta)^,A, = [Ta,A,], (15) 

where Ta and Ta are the d— and D— dimensional representations of the Lie algebra of the symme- 
try. Equations (fT4l) and ([T?] ) will be our guiding lines in defining states with prescribed symmetries. 



6. The Hamiltonian 

Given a matrix product state, the reduced density matrix of k consecutive sites is given by 

_ tr{{Al...Al®A,,...A,jE^-') 

The null space of this reduced density matrix includes the solutions of the following system of 
equations 

d-l 
Ji,--- Jfc=0 

Given that the matrices Ai are of size D x D, there are D"^ equations with d'^ unknowns. Since 
there can be at most independent equations, there are at least d'^ — solutions for this system 
of equations. Thus for the density matrix of k sites to have a null space it is sufficient that the 
following inequality holds 

d^ > D^. (18) 

Let the null space of the reduced density matrix be spanned by the orthogonal vectors \ea) , (a = 
1, ■ ■ ■ s, > — D"^). Then we can construct the local hamiltonian acting on k consecutive sites as 

s 

:= J^/ic.|e„)(e„|, (19) 

a=l 

where are positive constants. These parameters together with the parameters of the vectors |ej) 
inherited from those of the original matrices A^, determine the total number of coupling constants 
of the Hamiltonian. If we call the embedding of this local Hamiltonian into the sites / to / + A; by 
hi^i+k then the full Hamiltonian on the chain is written as 

N 



H = J2 Ki+k- (20) 
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FIG. 1 : (Color Online) A spin 1/2 ladder, an operator cr^^^ refers to the Pauli operator aa on rung k and leg 
number (i). 

The state | \E') is then a ground state of this hamiltonian with vanishing energy. The reason is as 
follows: 

AT 

{^\H\^) = tr{H\^){^\) = Y,HKi+kPi,i+k) = 0, (21) 

1=1 

where pi^k+i is the reduced density matrix of sites I iol + k and in the last line we have used the 
fact that h is constructed from the null eigenvectors of p for k consecutive sites. Given that is a 
positive operator, this proves the assertion. 

III. LADDERS OF SPIN ONE-HALF PARTICLES 

We now specify the above generalities to a ladder of spin 1/2 particles. The ladder consists of 
N rungs and obeys periodic boundary conditions. To each rung of the ladder we associate four 
matrices Aqo, ^oi; Aio and An respectively pertaining to the local states |0, 0), |0,1), |l,0)and 
|1, 1). Here we are using the qubit notation which corresponds to the spin notation as |0) = |+) 
and |1) = |— ). The first and the second indices refer respectively to the states of sites on the legs 
1 and 2 as shown in figure (1). The rungs are numbered from 1 to N with a periodic boundary 
condition. The operator cr^^ refers to the Pauli operator a a on leg number i and the rung number 
k. The total spin operator on a rung at site k is denoted by := + <y'k^)- 

A. MPS operators for observables 

Let us label the legs of the ladder by 1 and 2 as in figure (1). In each Aij the first and the second 
indices refer respectively to legs number 1 and 2. Any operator corresponding to leg number i is 
designated with a superscript [i). With these conventions, one can easily use ^ and write down 
the MPS operator corresponding to an observable. For example, for the magnetization in the x and 
z directions in rung 1 , we have respectively 

Em = Aoo ® Aio + Aio ® Aoo + Aqi ® An + An ® Aqi, (22) 

and 

Em = Aoo ® ^00 + ^1 ® ^1 - ^10 ® ^10 - An ® An- (23) 

z 
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For the total magnetization in the z direction in a rung we have 



Es^ = ^{E a) + E .2)) = {Aoo ® Aoo - Au ®An). (24) 
Other MPS operators can be obtained in a similar way. 

B. Symmetries of spin ladders 

In this paper we restrict the dimensions of our matrices to = 2, and demand that our models 
have time reversal symmetry so that the matrices Aij are chosen to be real. We will be looking 
for models which have rotational symmetry in the x — y plane of spin space. Thus we require that 
there be a matrix 7^ such that 

[%, Aoo] = Aoo, [%, An] = -Au, [%, Aoi] = [%, A^o] = 0. (25) 
It is an easy exercise to show that the solution of these equations is 

-"-(o!)' ^-(-)' ^-(n)' 

where T is found to be 7^ = | ^ ^ ^^'^ have excluded the solutions with g = or h = 0, 

which lead to trivial uncorrelated states. Hereafter we use the freedom in re-scaling the matrices 
(without changing the matrix product state) to set the parameter h = 1. Let us now consider extra 
discrete symmetries in addition to the above continuous symmetry. These are as follows: 

a: Spin flip symmetry represented by a matrix X, such that 

XA^^X-' = eAjj, (27) 

where e = ±1 and i means i + 1 mod 2. This symmetry imposes the following condition on the 
solution (|26|) 

a' = eb, b' = ea, (28) 

g' 

e 

b: Symmetry under the exchange of legs of the ladder represented by a matrix Y such that 



where X is found to be X 



YA,,Y-' = 7]A,,i, (29) 
where 77 = ±1. This symmetry imposes the following condition on (l26l) 

a = rja, b' = r]b, (30) 



where Y is found to be Y 



1 

7] 
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c: Parity or a left-right symmetry, represented by a matrix 11, such that 



liAi.U-^ = aAi,, (31) 



with a = ±1. 

This symmetry imposes the condition 



era, 



b' = aa, (32) 



on (|26|) where 11 is found to be 11 



1 
a 

Thus when we impose any of these discrete symmetries we are left with two three-parameter 

families of models, each family being distinguished by the discrete parameter e, r] or a. 

It is now readily seen that imposing any two of these symmetries makes the model symmetric 

under the third one too. A model which has all three symmetries is defined by the following set of 

matrices: 

Such matrices satisfy (|29l) with i] = ea. Thus equation [33] defines four two-parameter families of 
models on spin ladders which have 5*0(2) symmetry in addition to three types of Z2 symmetries. 
The families are distinguished by the pair of discrete parameters (e, a). 

d: Full rotational symmetry Let us now see if we can construct models which have full 
SU (2) symmetry, symmetry under rotations in the spin space. To this order we note that to have 
full rotational symmetry, the matrices defined by 

:= Aoo, 5i,o := ^(^01 + ^10), := A^, (34) 

and 

5o,o := ^(Aoi-Aio), (35) 

should respectively transform like the spin 1 and spin representations of the su{2) algebra, that 
is, we should have 

[Ta, So,o] = 0, a = X, y, z, (36) 

and 



[%, Bi^ni] = m Bi^rn, [Xt, 5/,m] = v2 - m(m ± l)fi/.m±i, / = 1, m = -1, 0, 1, 

(37) 

where Ta = ^aa, a = x,y, z form the two dimensional representation of the su{2) algebra. It is 
well known that the matrices ((T_, -^cTz, —cr+) transform like a vector under the adjoint action of 
SU{2). The matrix 5o,o should also be a multiple of identity. Thus we should set g = —1, and 
satisfy the following equations 

Aoi + Aio = cr„ Aoi - Aw = ul. (38) 
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This puts the constraints 

a + a = 1, b + b' = —1, a — a' = u, b — b' = u, (39) 
with the unique solution 

a = , a = , & = , b = — - — , (40) 

where u is an arbitrary real parameter. This will give us a one-parameter family of models with 
full rotational symmetry and spin-flip symmetry 



Remark: Note that comparison of these parameters with the constraints (l28l [30l) and (??) 
shows that full rotational symmetry is compatible with spin-flip symmetry for arbitrary values of 
the parameter u and compatible with parity or leg-exchange symmetries only for u = 0. 



The model (|4T1) has already been studied in [|l7h . To see the correspondence with that work, we 



can collect the above matrices in a vector-valued matrix A defined as .4 = Aij\i,j). In view 
of our notation for spins (|0) = |+), |1) = |— )) and the notation of [|l7n in which the single and 
the triplet states are denoted respectively by \s) and |to), the matrix A derived from 

(|4TI) becomes 

^^fu\s) + \to) -V2\h) \ 

v^V v^l^-i) u\s)-\t,)J'^ ^^^^ 
which modulo an overall constant is identical to the matrix given in fil7i1 . 



IV. PROPERTIES OF THE STATES AND CORRELATION FUNCTIONS 



We now study the properties of the states constructed above. For simplicity we consider in 
detail only two general classes. The first class is defined by (|33l) and has 5*0(2) symmetry in 
addition to three Z2 symmetries, and the second class is defined by (|4T1) which has full rotational 
symmetry in addition to one Z2 symmetry, the spin flip symmetry. As mentioned in the previous 
section, full rotational symmetry is compatible only with spin flip symmetry for generic values of 
the parameter u and is compatible with the other two symmetries only when u = 0. In order to 
derive results which can be specialized to the two classes mentioned above we study in detail the 
properties of the state, defined by the equation (l28l) . The matrices are now given by 



A 



01 



a 
b 



A 



10 



eb 
ea 



A 



00 



The matrix E for this class has the following form 

fa^ + b^ 

E = 



V 






1 





2ab 
2ab 




a 



9 




9 



2 + 



A 



11 





1 



(43) 



(44) 
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with eigenvalues 



Ai = + 6^ + = + 6^ - g, A3 = A4 = 2ah. 

For > 0, the largest eigenvalue is Ai and for g < it is A2. Hence the point (7 = is a point of 
phase transition. The right and left eigenvectors of E are simply obtained and one can determine 
all the relevant quantities of the ground state in closed form, in straightforward way after some 
rather lengthy calculations. The reduced one-rung density matrix is obtained from (??): 











\ 



0^ + 52 2eah 
2ea6 ^ 6^ g 
\0 ^ \g\J 



(45) 



where Q := 2d^ + 2}? ^2\g\. This matrix can be rewritten as 



P = Q O^l(|00)(00| + + 2ea6(|01)(10| + |10)(01|) + (a^ + 6^)(|01)(01| + |10)(10|)] , 

(46) 

(47) 



or 



where we have used the notation introduced in equation P^ . 



From this density matrix one can obtain a lot of information about the observables pertaining 
to a single rung. For example it is readily seen that the average magnetization at each single site 
and hence each single rung is zero, i.e. 



(2)\ 



(S) = 0, 



(48) 



where S is the spin of a rung, implying an anti-ferromagnetic state in which every single site is in 
a completely mixed states. It is also seen that 



\g\-a^-h'' 
\g\ + a? + h^'- 



2eah 



(49) 



where n is any unit vector in the x — y plane. 

Defining the total spin of a single rung as S := + ff^"^^), we find from the above result that 



{a^ehf + 2\g\ 
a? + b^ + \g\ ■ 



(50) 



Thus for \g\ 00, each rung will be in a mixture of spin one states, but for g = 0, the spin one 
and spin zero multiplets can mix, depending on the value of a and b and e. The entanglement of 
this rung with the rest of the lattice is measured by the von-Neumann entropy of this state, defined 

as S* = —tr(p log p). 
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From the eigenvalues of the one-rung density matrix, 



\9\ {a - bf (g + bf 

ai = a. = ^, «3 = ^^, "4 = ^^. (51) 



we readily find 



S = \ogQ-^[\g\ \ogg^ + (a + hf log(a + hf + (a - hf log(a - hf] . (52) 

As ^ oo 5 ^ log 2 = 1. Therefore in this limit each rung is still entangled with the rest of the 
ladder. This means that a rung is not in a pure state and the state of the ladder is not a product of 
single rung states. 

The entanglement of the two spins of a single rung with each other is measured by the concurrence 
^ of this state which for the density matrix (|45l) is given by 

C = max{0, 2amax - 1), (53) 

where ctmax is the largest eigenvalue of the matrix p. A careful analysis of the eigenvalues shows 
that 

2|a6| - 1^1 
a2 + 62 + 

As ^ oo, C ^ which means that although each single rung is not a pure state, it is a 
separable state. In fact these can also be verified directly by looking at the one-rung density matrix 
in this limit. From (|46|) we have 



C = max{0, ^\ ,'^^, ). (54) 



P(^7^oo) = ^(|00)(00| + |ll)(ll|). (55) 

Finally we calculate the correlation functions of different components of spins of the rung as a 
function of the distance between the rungs. We find from ^ that 

{S.,S.,.) = -9'^-^^^^^ = {S.A,)e~^^-^y^^, (56) 



with a longitudinal correlation length 



• , , „2i^2_i_i„i , ) (57) 



and 



(^n,i^n,.) = {sgn{g) + e)(a + sgn{g)hf^--4^^§^^^^ = {S^,,S^,2)e-^^-'y^- (58) 



2 (a^ + b^+\g\) 

with a transverse correlation length 



1 / g^+fe^+lgK - ^^^^ 
|2afe| ) 
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— 1^ 
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FIG. 2: (Color Online) The longitudinal and transverse correlation lengths of matrix product states |43J as a 



function of the parameter x for different values of the parameter /i := 
and longitudinal correlation lengths diverge at x = 0. 



2|afc| 



At ^ = 1 both the transverse 



It is seen that the longitudinal correlation length depends on a single parameter x := as 
= In(i^q^) and the transverse correlation length depends on x and another parameter fi : = 

l^p^ < 1 as = ln(ii|^). Figure (2) shows the behavior of these correlation functions for 
different values of the parameters x and yU. 

An interesting point about transverse correlation functions is that depending on the sign of e, 
it becomes identically zero, on one side of the g— axis and different from zero on the other side. 
Thus if one insists on a definition of an order parameter to be zero in one phase and non-zero in 
the other, then we can safely say that in these models, the transverse correlation function is an 
order parameter which signals a quantum phase transition. 



Let us see study some limiting cases in these general models. At the point of phase transition 
5^ = 0, as seen from (|46|) . the state of a rung becomes a mixture of spin-zero states, and the two 
spins of a rung become fully anti-correlated in the z direction, a fact which is reflected in (|49l ). Near 
this point the correlation length becomes very large as seen from (l57l) . although the amplitude 
becomes small since it is proportional to \g\. Thus at g = there is no long-range order in the 
model. Also from (|54l) and (ISTI) it is readily seen that at g = 0, the concurrence (or entanglement) 
of the two spins of a rung becomes maximum and equal to 



Cig = 0) 



2\ab\ 



(60) 



On the other hand from (l52l) we see that as we approach the point g = from both sides, the 
von-Neumann entropy decreases. In addition, the derivative of both types of entanglement become 
singular at g = 0. These facts show that the point of phase transition in these systems, is a point 
where the spins of a rung become highly entangled with each other and each rung becomes only 
slightly entangled to the rest of the lattice. 

We can also obtain the explicit form of the state in this limit. From (|43l) . and ^ we see that in 
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this limit, no two spins in a rung can be in a ^ |J ^ state. A little reflection shows that they can not 

be in the state ^ | ^ either, since in the absence of Aqq, the only string of matrices with non-zero 

trace is a string of matrices Aqi and Aiq in arbitrary order. Since these matrices commute with 
each other, the resulting state has a simple description. Define two local states of a rung as 

where for convenience we have re-introduced the spin notations + and — instead of qubit notation 
and 1. Define a global un-normalized state d^~^) to be the equally weighted linear combi- 
nation of all states with k local u states and N — k local d state. Then from (|43l) and (O we find 
that the ground state of the chain in the limit g = Ois given by 

1^) = [a'^(e6)^-'= + 6'=(ea)^-'=] (62) 

where Z is the normalization constant given by 

Z = J2[ k) [a'i^bf-'^ + b'iea)^^']' = 2 [{a^ + bY + i^abf] ■ (63) 

k=0 ^ ^ 

At the other extreme when \g\ — > oo, the state of each rung becomes a mixture of fully aligned 
spins either in the positive or negative z direction. This is also reflected in (|49l) . In this limit a rung 
becomes entangled with the rest of the lattice, sine lim|g|_,oo S{p) = 1 and the two spins of a rung 
become disentangled from each other since C = 0. 

The explicit form of the state can also be obtained in this limit. In this case the only strings of 
matrices with non-vanishing traces are strings of Aqq and An in alternating order. Thus if we 
define two local states 

then the ground state in this limit will be a GHZ state 

In the next subsection we specialize these results to the two classes we discussed in the beginning 
of this subsection. 



A. Class A 

This is the class which has the S0(2) symmetry (rotation around the z axis in spin space) and 
three Z2 symmetries (spin flip, parity, leg exchange). For this class we have from (|32|) that b = era. 
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FIG. 3: (Color Online) The von Neumann entropy of one single rung and its concurrence in models A, as 
functions of the parameter x := which is taken to be dimensionless. 



The parameter x defined after equation (1591) becomes equal to 

. := (66) 

Insertion of 6 = era in various quantities of the previous subsection shows that all the quantities 
can be expressed as a function of x, namely we find 



and 

{Sn,iSn,r) = {l + esgn{x)){l + a sgn{x))\x\-ij--^j^, ^" " h^^TT^' 



We also find 

5 = -^^(l-log|x|) + log(l + |a;|), (70) 





X 




1 + 


x\ 



and 



C = max{0,- (71) 



1 + \x 

Figure (2) show the entropy and the concurrence of the state of a single rung for this model. 



B. Class B 



This is the class which has full S0(3) symmetry (rotation in spin space) and one discrete sym- 
metry (spin flip) for generic values of u and all the three Z2 symmetries for m = 0. For this class 
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FIG. 4: (Color Online)The von Neumann entropy of one single rung and its concurrence in model B as 
functions of the parameter u, which is taken to be dimensionless. 



we have from (|4TI) g = e = —1, a = and b = Inserting these values in the equations of 
the previous section we find the following: 



and 



= -— (72) 

6 + 



(5n,l5n,.) = -4 \ „ , (73) 



where n is any direction in the spin space. The eigenvalues of the one-rung density matrix (1511) in 
this case will be 

1 

+ 3 M"' + 3 

Thus we find 

S = log(«^ + 3)-^^^, (75) 

and 

- 3 

C = max(0,^— ). (76) 
+ 3 

Figure (3) show the entropy and the concurrence of the state of a single rung for this model. 



V. DISCUSSION 

Although ladders of spin 3/2 models has already been studied in the context of vertex state 
models in lll2ll . in which the ground state is constructed by a suitable concatenation of vertices as- 
signed to single sites, the method developed in this work, in which each single rung is considered 
as a single site in a hyper-chain and the ground state is constructed as a matrix product state seems 
to us as more powerful and very easy to generalize to other spin models. In this work we have 
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applied the matrix product formalism for construction of models on spin 1/2 ladders. These mod- 
els have been constructed to have special discrete and or continuous symmetries and to display a 
quantum phase transition in a broad sense, that is displaying non-analytical behavior in their corre- 
lation functions. Naturally these non-analytical behavior can also be observed in the entanglement 
properties of pairs of spins in these models. This route can be followed to develop other models, 
e.g models with higher spins or alternating spins on the rungs, with alternating coupling constants, 
or next-nearest neighbor interactions. By exploiting higher dimensional matrices and having more 
free parameters at our disposal, we may be able to construct continuous families of models on 
spin ladders which have full rotational symmetry. In this article we have restricted ourselves to 
frustration-free Hamiltonians and MPS states with fixed-size. Relaxing this latter condition allows 
one to represent any state G C"'^^ as a matrix product state 10, Ea] and then one may be able 
to study more diverse kinds of phase transitions on spin ladders. These matters will be taken up in 
separate publications. 
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VII. APPENDIX 

In this appendix we briefly discuss the derivation of the explicit form of the Hamiltonian in 
terms of local spin operators. For simplicity we consider models in class A. The Hamiltonian for 
models in class B, can be constructed along similar lines. As explained in the main text, the starting 
point is to solve the system of equations 

1 

CijkiAijAki = 0. (77) 

i,j,k,l=0 

This solution space is 12 dimensional since we have 4 equations for 16 unknowns. Thus we should 
find a set of 12 orthogonal vectors which span this solution space and then form a non-negative 
linear combination of the corresponding one dimensional projectors. The Hamiltonian constructed 
in this way does not necessarily have the symmetries imposed on the state, unless we choose 
new linear combinations of these vectors which transform suitably under the symmetry operators. 
These new basis vectors are found to be: 

|2,2) := |0000) 

|2, 1) := i [-e|0001) + |0010) -e|0100) + |1000)] 

|2, 0) := ^ [-40^(10011) + |1100)) + g (|0101) + (Te|0110) + cre\1001) + |1010))1 

|1, 1) := i [e|0001) + |0010) -e|0100) - |1000)] 
|1,0) := ^ [|0110) - |1001)] 
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1) := ^ [-(xelOOOl) + |0010) + e|0100) -(t|1000)] 
|1',0) := -i=[|0101)- 11010)] 

|0,0) := ^ [|0101) -(xelOllO) -ae|1001) + |1010)], (78) 
with |/, — m) = erf |/, m), V l,m. 

Here we have organized the states in multiplets, with a labeling reminiscent of the one used in 
labeling the states of so(3) representations. The reason is that in a certain limit (e = a = g = 
— 1, a = |) these states actually comprise the irreducible representations of so(3) which emerge 
from the decomposition of the product of four spin one-half representations living on the sites of 
two rungs of the ladder. In fact for the product of these spin 1/2 representations we have 




= 2 © 1 © 1' © 1" © © 0'. 



However only the representations 2, 1, 1' and span the solution space of (l77l) in this limit. This 
labeling is useful because we can see under what conditions, the Hamiltonian become fully rota- 
tional invariant. Note also that we have used the labeling \i j k I) in accordance with the labeling 
of (1771) for Cijki (see figure (5). The reader can easily check that each of the above states (or 
more precisely the corresponding one dimensional projector) is invariant under the parity opera- 
tion 11: \i j k I) \k I i j) and leg-exchange of the ladder Y : \i j k I) ^ \j i I k) . They also 
transform to each other under the spin-flip operation X : \i j k I) \i j kl) . 

The local Hamiltonian which has the three discrete symmetries mentioned in the text and the 
symmetry around rotation along the z axis in spin space, is constructed as follows: 

I 

h= ^ '^nim{\l,m){l,m\ + \l,-m){l,-m\), (79) 

Z=2,l,l',0 m=0 

where the coefficients f^i^m are non-negative parameters and together with the parameters a and g 
form the 10 free parameters of the Hamiltonian. Of course the total number of coupling constants 
(interaction strength) is 8, since we can always shift the ground state energy and also set the scale 
of energy by redefintion of these parameters. Re-expressing the local operators in terms of Pauli 
spin operators and rearranging terms, we find, after a rather lengthy calculation, the Hamiltonian 
acting on the ladder where we have multiplied hhy a factor of 8 for convenience, (see figure (6) 
for labeling of sites of the ladder). Note that we use zt as an abbreviation to denote the operator 
ai^z in the following. The total Hamiltonian is then given by 

H = Hi + H2 + + + H['f (80) 

where 

Hi = yjo+Ji 
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FIG. 6: (Color Online) The labeling of sites of the ladder for reading the interaction terms in the Hamilto- 
nian. 

H2 = ^ Jo + Ji ^i'^i'+i + J2 o-i' ■ (81) 

i 

i 

Hi2 = -^3 0"j ■ + J4 (O^i • + ■ + J5 ■ ■ ^i'+l) 

i 

+ Jq {ffi ■ ai+i){ai' ■ ai'+i) + J7{ai ■ ai'+i){ai' ■ ai+i), 

and 



^12 = "^11 {ziZi>{ai+i ■ ai'+i) + Zi+iZi>+i{ai ■ ai')) 

i 

+ J12 ■ O-j'+l) + ■ fj+l)) 

+ Ji3 {ziZi>+i{ai ■ ffi+i) + Zi'Zi+i{ai ■ ffi'+i)) . 

Thus there are bond interaction and plaquette interaction in the Hamiltonian all depending on 10 
parameters. The coupling constants are related to these parameters as follows: 



Jo = 1^22 + 4(/i2i + /in + /ii'i) + 1^10 - l^i'o + 2(/ioo + 1^20) + 16aV2o 
Ji = 1^22 + 7:(-At2i + Atii + o-jiyi + /ii/Q - /iio) + -a^igae - 2a?)^2o 



2 

J2 = -{fJ'21 - AHi - <7fJ'i',i> - -o-ea gfi2o 

,2 



J3 = -e(/i2i - /ill) - o-e/ii/i + eo-(^5fV20 - /^oo) 
J4 = --(^{fJ'21 + /ill - fJ'i'i) - -a 91^20 
J5 = ^(/ioo - /iio - /ii'o) + ^/i2o(5'^ - Sa"^) 

Jg = 20^/^20 + 2 (^I'o ~ /^lo) 

8 1 
J7 = :^o,'^l^20 + ^{fJ'io - f^i'o) 

2 

Jg = 2/^22 + e(/i2i - /ill) + cre/ii'i - /iio - l^i'o + (c^e - l)/ioo + 2/^20(80'' - (1 + ae)/) 
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^9 = A*22 + ^e(/X2i + //ii - i^i'i) + ^(/Uio - fJ'i'o) + ^a^(5' - 2a^)//2o 

2 8 

^10 = A^22 + 2(e - l)/i2i + (e + 1)((7 - l)^i/i + (1 - (7e)^oo + ^(1 + cre)g^iX2Q + -a^(2a^ - - gcre)iJ,2o 
Jii = ^(-eAi2i + e/xii - ea/^i/i + /^lo + /xi'o) + ^(ecr - l)/xoo + ^(Sa"^ - (1 + o-e)5'^) 
Ji2 = ^(a«2i - A«ii - cr/iri + /iio - /ivo) - ^a^(2a^ - gae)ii2o 
Ji3 = - e/iii + e/ii/i - /xio + firo) + ^0^20(5' - 2a^). 

This is a complicated looking Hamiltonian, with many types of interactions, but in view of the 
large number of parameters, it is possible to look at specific subsets of the parameter space, where 
some of the interactions are absent. Indeed the parameter space consists of four disconnected 
parts, each of which corresponds to one choice of the pair (cr, e). Let us for example consider the 
subset on which the Hamiltonian has full rotational symmetry. It is well known that any operator 
of the form Ylm=-i 1^' "^)(^' "^1 ^ scalar. Thus in the limit g — e — a — —1, a — | if we set 
A<2m 6//, ni^m —'■ 2i^, Hv,m —'■ '^^ V m and //o,o 277, then we will have a Hamiltonian which 
is fully rotational invariant. On can see that in this limit, all the coupling constant corresponding 
to non-scalar terms in the Hamiltonian vanish and we are left with 

Jo = 48/x + lOi/ + 6^ + 4// 

J2 = 5// - 1/ + ^ 

J3 = 10/x - 2z/ - 2^ - 2?7 

J4 — 5fi + u — C, 

J5 = fl-U-^ + T] 

Jq = + ^ - z/ 
J7 = ji + v-S,. 
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